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ABSTRACT 

The spectra obtained above 100 MeV by the EGRET experiment aboard the Compton Gamma-Ray 
Observatory for a handful of gamma-ray bursts has given no indication of any spectral attenuation that 
might preclude detection of bursts at higher energies. With the discovery of optical afterglows and 
counterparts to bursts in the last few years, enabling the determination of significant redshifts for these 
sources, it is anticipated that profound spectral attenuation will arise in the Gamma-Ray Large Area 
Space Telescope (GLAST) energy band of 30 MeV-300 GeV for many if not most bursts. This paper 
explores time-dependent expectations for burst spectral properties in the EGRET/GLAST band, focus- 
ing on how attenuation of photons by pair creation internal to the source generates distinctive spectral 
signatures. The energy of spectral breaks and the associated spectral indices provide valuable informa- 
tion that constrains the bulk Lorentz factor of the GRB outflow at a given time. Moreover, the distinct 
temporal behavior that is present for internal attenuation is easily distinguished from extrinsic absorp- 
tion due to intervening cosmic background fields. These characteristics define palpable observational 
goals for both spaced-based hard gamma-ray experiments such as GLAST, and ground-based Cerenkov 
telescopes, and strongly impact the observability of bursts above 300 MeV. 

Subject headings: gamma-rays: bursts — radiation mechanisms: non-thermal — gamma rays: theory 
— relativity 



1. INTRODUCTION 



^ c~| High energy gamma-rays have been observed for six gamma- 
fffry bursts by the EGRET experiment on the Compton Gamma- 
lay Observatory (CGRO). Most conspicuous among these 
gipservations is the emission of an 18 GeV photon by the 
+^RB940217 burst (Hurley, et al. 1994). An additional case of 
interest is provided by the so-called MILAGRITO burst (Atkins 
;et al. 2000; Atkins et al. 2003), GRB 970417, with its uncor- 
. Jn)borated detection of ~ 3cr significance in the TeV band. Tak- 
S^g into account EGRET 's field of view, its detections indicate 
jujiat emission in the 1 MeV-10 GeV range is probably common 
Qimong bursts, if not universal (see Dingus 1995 for a discussion 
of EGRET high energy burst statistics). One implication of 
GRB observability at energies around or above 1 MeV is that, 
at these energies, spectral attenuation by two-photon pair pro- 
duction (77 — > e + e~ ) is absent in the source. From this fact, 
Schmidt (1978) deduced that a typical burst had to be closer 
than a few kpc, if it produced quasi-isotropic radiation. 

In the aftermath of BATSE's revelation (e.g. Meegan et 
al. 1996) that most if not all long bursts are at cosmological 
distances, Krolik & Pier (1991) and Fenimore et al. (1992) 
proposed that GRB photon angular distributions are highly 
beamed, being produced by a relativistically-moving plasma, 
a suggestion that has become an underpinning of the GRB 
paradigm. This can dramatically reduce r 77 below the ~ 10 11 — 
10 12 values realized for isotropic photons, and blueshift spectral 
attenuation turnovers above the observed spectral range. De- 
terminations of the bulk Lorentz factor T of the GRB medium 
have mostly concentrated (e.g. Epstein 1985; Krolik & Pier 
1991; Baring 1993) on cases where the angular extent of the 
source was of the order of 1/F. These calculations generally 
assume an infinite power-law burst spectrum, and deduce (e.g. 



Woods & Loeb 1995; Baring & Harding 1997b) that gamma-ray 
transparency up to the maximum energy detected by EGRET 
requires r ;> 100 -10 3 for cosmological bursts. Since BeppoSax 
spawned the age of precise determination of GRB redshifts, such 
bounds can now be refined for the BeppoSax, HETE and Swift 
databases; the discovery of high redshift bursts (e.g. see Cum- 
mings et al. 2005; Kawai et al. 2005; Berger et al. 2006; Haislip 
et al. 2006, for Swift bursts) might suggest even higher esti- 
mates of bulk Lorentz factors in those sources. 

While the power-law source spectrum assumption is expedi- 
ent, the spectral curvature seen in virtually all GRBs by BATSE 
(Band et al. 1993; Preece et al. 2000) is expected to play an 
important role in reducing the opacity for super-GeV and TeV 
band emission from these sources (Baring & Harding 1997a). 
Such curvature is patently evident in 200 keV-2 MeV spectra 
of some EGRET-detected bursts (e.g. Schaefer, et al. 1992), 
and its prevalence is indicated by the generally steep EGRET 
spectra for bursts (e.g. Hurley, et al. 1994; Schneid, et al. 
1992; Sommer et al. 1994). This curvature profoundly im- 
pacts 77 opacity determinations, yielding a huge dearth of tar- 
get soft photons in comparison with power-law spectra extrapo- 
lated down to the classic X-ray band. Hence, when considering 
emergent photons of dimensionless energies e 7 > T 2 /e B , where 

e B is the dimensionless energy (in units of m e c 2 ) of the spectral 
break in the BATSE band, the target photon energy is below 
e B m e c 2 so that accurate pair opacity calculations mandate de- 
tailed treatment of the spectral curvature observed in bursts. 

This paper enunciates the principal properties of pair 
production opacity that couple to spectral shape in the 
BATSE/EGRET energy range, embellishing upon and focusing 
the work of Baring & Harding (1997a), so as to identify possible 
observational diagnostics for the Gamma- Ray Large Area Space 
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Telescope (GLAST; http://glast.gsfc.nasa.gov/) mission 
and ground-based experiments such as the MILAGRO, MAGIC, 
HESS and VERITAS atmospheric Cerenkov telescopes. More- 
over, it helps define GRB science drivers for future initiatives in 
the hard gamma-ray band above 100 MeV. These spectral signa- 
tures at energies where significant opacity is realized are clearly 
distinguishable from absorption by background radiation fields, 
and in particular via their time-dependent evolution. The fact 
that turnovers produced by constant Lorentz factor evolution 
can be easily discriminated from those generated by decelerat- 
ing expansions is compelling. It may also be possible to elicit 
details of the adiabaticity or otherwise of the post-fireball ex- 
pansion. Accordingly, hard 7 -ray telescopes can, in principal, 
provide powerful probes on such evolution of bulk motions in 
bright bursts, which can then be used to generate useful bounds 
on the explosion energy, a key piece of information elucidating 
the nature of the central engine for bursts. Before discussing 
these evolutionary effects and the potential for observational di- 
agnostics in Sections 13 and 01 the paper presents refinements in 
the formalism for pair production opacities in Section [21 

2. GAMMA-GAMMA PAIR PRODUCTION OPACITIES 

In this Section, a detailed formalism for 77 pair creation 
opacities from piecewise-continuous power-law distributions of 
photons is expounded upon, extending the presentation of Bar- 
ing & Harding (1997a) by providing greater detail and also com- 
pact analytic reduction of the integrals for the optical depths. 

2.1. General Formalism for Internal Target Photons 

To assess the role of two-photon pair production in attenu- 
ating gamma-ray burst spectra in the EGRET/GLAST band, 
the interactions of photons created only within the emission re- 
gion are considered here, neglecting the presence of any external 
radiation; the impact of such external radiation fields on burst 
spectra is discussed in Section l2~4*l below. The calculations of this 
subsection mirror those of Baring (1994), whose work is used as 
a basis for elucidating the effects of GRB spectral structure on 
pair production opacities. For simplicity, the photon distribu- 
tion is assumed to be azimuthally symmetric about some beam- 
ing axis (that is presumably aligned near to the line of sight to a 
distant observer), with the spectrum n(e) being identical for all 
photon angles 9 £ relative to the axis within some cone of emis- 
sion. Here e is the photon energy in units of m e c 2 , a convention 
adopted throughout this paper. The photon distribution func- 
tion then takes the form n 7 (e, fi s ) = n(e) f(p £ ) , with f(fJ> e ) 
normalized to unity, and where fi £ — cos 9 £ is the cosine of the 
photon angle with respect to the axis of symmetry. Further, for 
anisotropic (beamed) photons we take the simple case that the 
angular distribution about the axis is uniform within a cone of 
half-angle 9 £ = 9 m : 

/(Me) = i _ 1 „ , Mm = C0S6> TO < H e < 1 . (1) 

r Mm 

This conical sector form is not exactly equivalent to the angular 
distribution of an isotropic photon distribution boosted by some 
bulk Lorentz factor T ~ l/0 m with respect to an observer at 
infinity, since it lacks the "wings" at fi £ <; Mm . Yet it suffices 
to exhibit the general nature of photon beaming effects on pair 
production rates. Consequently, the results presented here, for 
example that inferred from Eq. (|17fl . differ slightly from what 
would be derived from a beamed adaptation of the isotropic re- 
sults of Gould and Schreder (1967, as expounded briefly below), 



which was the approach of Krolik & Pier (1991) and Fenimore, 
Epstein & Ho (1992). Scientific conclusions pertaining to rel- 
ativistic beaming in GRBs are not contingent upon such sub- 
tleties, largely because the computed optical depths are most 
sensitive to the beaming factor 9 m , i.e. 1/r, among the rele- 
vant parameters. Moreover, adherence to one particular angular 
distribution is not mandated by observations nor by theoretical 
insights. For example, diffusive particle acceleration at the rela- 
tivists shocks that are suggested as being sites for energization 
in gamma-ray bursts does not generate isotropic distributions 
at the shock in any frame of reference (e.g. see Bednarz & Os- 
trowski 1998; Kirk et al. 2000). 

For photons that interact with themselves to produce pairs, 
a slight modification of the form for the pair production op- 
tical depth r 77 (e) in Eq. (7) of Stepney and Guilbert (1983) 
was adopted by Baring (1994) as a starting point for analytic 
developments. This form, appropriate for photon distributions 
independent of azimuthal angles, is 

7Vf( e ) = J d^efiPe) J dfj, u f(n u ) J duj^^- 

(2) 

x j X+ X 3 o-^(x)d X 

•'x- \J (x 2 - x 2 ) (x 2 - X 2 ) 

Here \ = — cos0)/2] 1 / 2 is the center-of-momentum (CM) 
frame energy of the photons, whose energies and angle cosines 
in the observer's frame are (e, Me) and (u>, /i w ) respectively. 
is the angle between the photon directions. The pair pro- 
duction threshold condition is then x > 1 . Since \x £ = cos 9 £ 
and fj, u — cos 9 U , the range [0, n] of azimuthal angles in the 
observer's frame establishes the bounds 

X 2 = ^{l-cos(0 E ±*y} . (3) 

The focus of Baring (1994) was on infinite power-law number 
density distributions n(uj) . In this paper, low energy cutoffs 
with u) > will be introduced, since these prove expedient 
in treating piecewise continuous broken power-laws, which con- 
stitute suitable approximations to the general shape of burst 
spectra. Specifically, the form 

/ s /0 , e < e c 

is used for some minimum photon energy e c m e c 2 , following 
Gould and Schreder (1967) and Baring & Harding (1997a); note 
that the works of Baring (1993, 1994) used an integral photon 
flux index. The units of both n 7 and n(e) correspond to inverse 
volumes. The change of variables ( = (1 — cos0)/2 = x 2 ' jeu) 
can then be employed to facilitate the manipulation of the inte- 
grals along the lines of Baring (1994). Defining a parameter 

v = i ( 5 ) 

that is the center-of-momentum frame energy for head-on pho- 
ton collisions, the form 

r 77 (e) = n^RTa^Om , ^feeZ^j e a ~ x , (6) 

is obtained, where 

T a (9 m , r)) = ± r dx^M- F a (0 m , V, X) (7) 
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is an integration over the pair production cross-section <r 77 , 
the angular distributions contribute the factor 



Fa{Qm, V, X) 



1 



n J» m 1 - V™ 7 Mm 1 - Mr, 



and 



(8) 



c- V(C+-0(C-C-) Vx 

Here, Q(x) is a Heaviside step function such that Q(x) = 1 
for < x < 1 and is zero otherwise. It expresses the kine- 
matic condition \ 2 > ££cC f° r C — (1 — cos9)/2 that is im- 
posed by the spectral truncation. The £ -integration limits are 
C± = X±/ £UJ = [1 — cos(9 E ± 9 w )]/2 . Note that in the case of 
e c = , effecting the substitution a — > a + 1 in these results 
reproduces those of Baring (1994). 

It is elucidating to focus first on the case of isotropic photons, 
i.e. when 9 m = n . The angular integrations in Equation JSJl arc 
facilitated by following Baring (1994) and changing variables to 
with 2£ = {C + + £- + (C+ ~~ C-) cos0} , and then perform- 
ing the solid angle transformation 4>) — > (p 6 u>, 4>eu) , where 
P"eu> = cos0 is the cosine of the angle between the photons. 
This results in a dramatic simplification of the three integrals. 
Defining 



min< 1 



X 
V 

Equation JSJ becomes, in the limit 9 n 



V, x) = 



q 



2(a+l) 



1 



(9) 



(10) 



with the rj < 1 case (i.e. infinite power-law; q = 1 ) reproducing 
Eq. (10) of Baring (1994). It then follows that 



v) a + 1 



(11) 



2.2. Angular Contributions for Strong Beaming 

Exploration of the angular integrals for cases \x m ^ — 1 , and 
specifically those where beaming is strong and /_t OT w 1 , is not as 
straightforward as the isotropic photon situation. Equation 
must be manipulated considerably to extract compact analytic 
forms. In the limit 1 — \i m ss 9^/2 <C 1 , to leading order in 
(1— fi m ) , it is easily shown that C± ~ ( a/1 — [x s ±y/l — ) 2 /2 , 
suggesting the change of integration variables to k such that 
C = k 2 (1 — ji m )/2 together with 1 — jj, £ = (1 — /i m )x 2 and 
1 — fi w = (1 — pbmjX^ . This yields the following form for : 



^(^m, ??, X) 



i6 e 

2a+i ^ 



g dx g 

^0 



K _ v/[( K+ ) 2 -«2] [ K 2_( K _)2] 



e 



2X 



(15) 
« 1, 



for fc-t = x e ± x w . Here, the fact that the integration is sym- 
metric under the interchange x £ <-» x w has been used to write 
2?uj < a^e without loss of generality, thereby introducing an extra 
factor of two. 

Before evaluating Eq. (|15fl in generality, it is instructive to 
identify compact results in two asymptotic regimes. In the par- 
ticular case of infinite power-law spectra, where rj = , the value 
of T a is independent of \ so that the angular integrations sep- 
arate from the energy ones. Then, further changes of variables 
along the lines of x u = wx e and k 2 — x £ [(l — w) 2 + 4wz] leads 
to evaluation of the integrals using identities 3.197.3, 9.134.3 
and 7.512.4 from Gradshteyn & Rhyzik (1980): 



A(a) 



A(a)9% , 

2 1 ~ 2a r(2a + 2) 
r(a + 2)r(a + 3) 



(16) 



where 



H(a, rj) 



g2(a+l) 



'11 



(x) d X 



(12) 



This completely defines the pair production optical depth for 
the case of an isotropic distribution of photons with a power- 
law spectrum that is truncated at low energies, a result that is 
quickly shown to be equivalent to Eqs. (21)-(23) of Gould and 
Schreder (1967) by observing their notation s — ► \ 2 > an d by 
performing the appropriate integration by parts. 

An exact expression for T a {-n, 1) was obtained by Svens- 
son (1987; see his Equation [B6]), and can be approximated 
(Baring 1993, 1994) to better than 1% for 1 < a < 7 using 
7{(a, 1) R* 7/6/a 5 / 3 . An asymptotic form for 77—5-00 can 
be obtained most expediently by adopting a form equivalent to 
Equation (jTT|) : 



V = — 



X 



2a+3 



X 3 t7~ 



, (x) dx 



(13) 



Using the asymptotic form a 11 (x) w (3cr T /8) [21og e 2x — l]/x 2 
for x ^> 1 , i.e. for the extreme Klein-Nishina domain, one 
quickly arrives at 



3 11, _ 1 , 1 
log e 2n - 1 



2a 



T 



2a 



r) > 1 ■ (14) 



This is identical to Eq. (24a) of Gould and Schreder (1967). 



valid for r\9 m jx *C 1 • This result is precisely that in Eq. (12) 
of Baring (1994) after the adjustment for the different spectral 
index convention used there; Baring (1994) also posited an ap- 
proximation, A{a) w 2/ (4/3 + a) 27 / 11 that is accurate to better 
than 1% on l<a<7. Equation l|16|) implies 



T a {9 m , r)) Rj A(a)H(a, 1)9 



2<i 



r]9 r , 



< 1 



(17) 



This can be inserted into Eq. © to obtain the corresponding 
final result for the pair production optical depth for cases where 
the beaming is strong, but when the spectral cutoff of the pho- 
ton distribution does not come into play. 

The opposite asymptotic limit arises for i]9 m /x S> 1 , and is 
readily tractable since it samples k < 1 and x £ rj x^ domains. 
This results in 



F a (0 m , v, x) 



o2n 



1 



X 

r\9„. 



2(a+l) 



r}6„ 



X 



> 1. (18) 



indicates that the coefficients for 
3> X domains differ, largely due 



Comparison with Eq. (|16fl 
the n9 m <C x and r\9 m 
to the behavior of T a near the critical kinematic value of 
V@m = X ■ Observe that this asymptotic domain establishes 
Fa{9m, T), x) K 1/^m > which is essentially a phase space factor 
in the angular integrals. 



4 



To reduce Eq. (|15f) for general values of rjO m /x , a moderate 
amount of manipulation is required. An effective approach is to 
change the integration variables via x^ = wx E and k = x £ p , 
so that x £ , w and p constitute the new set of variables. The 
w integration is then performed first, leading to the appearance 
if a simple inverse trigonometric function. The procedure is 
detailed in the Appendix, where it is found that T a is express- 
ible in terms of elementary functions and the hypergeometric 
function F = 2F1 . The attractively compact result of these 
manipulations is (for 9 m <C 1 ) 



A(a) 



* > 1 



, < * < 1 , 
where = x/Mm) , and Q a (^) is given by Eq. (|X4|) : 



(19) 



&«(*) 



a ) j arccos* - * 2 



7T ] a + 1 



a + 2 



#F(l/2, a + 3/2; a + 5/2; W 2 ) 
(a + l)(2a + 3) (a + 2) 



(20) 



Accordingly, much of the character of T a is encapsulated in 
the one parameter \I/ . Numerical protocols for the evaluation 
of the hypergeometric function F(l/2, a + 3/2; a + 5/2; \f 2 ) 
are outlined in the Appendix. Observe that Ga(^) < A(a) for 
< $ < 1 . Using this formalism, Eq. Q can now be written 



T a {6m, v) 



3 2n 



/C(o!, 77^ 



(21) 



where 



/C(a, 77) 



<7 T 



Hv grr(x) 

"X 2a- 1 

A 



.4(a) 



J g 77(x) 
«X 2a-l 
A 



v} 



(22) 



The general character of /C(a, 77) is illustrated in Fig. ^ in- 
cluding its comparison with the similar function Ti.(a, 77) that 
pertains to the isotropic photon case above. 

Further attempts to reduce the integration over \ m the ex- 
pression for T a lead to undesirable algebraic complexity, intro- 
ducing somewhat obscure special functions that preclude further 
insight. Hence numerical evaluations of the integral over \ suf- 
fice. These can be checked using two identifiable asymptotic 
formulae. The first is the r)9 m < 1 form that is encapsulated 
in Eq. l(T7l) . which uses K.(a, 77) —> JC(a, 1) = A(a)H(a, 1) for 
i] < 1 . The second is for the extreme Klein-Nishina limiting case 
of rjO m ^> 1 , which can be obtained by inserting Eq. I|19fl into 
Eq. (JJJ , and then invoking the integral evaluations in Eq. (|A9|> . 
The result is 



1~ a (6 m , 77) 



i^{ l0 Se^m-! + ^}, 

(23) 

which provides an excellent approximation to l~ a (9 m , 77) for 
V@m <; 4 . The similarity of this to the Klein-Nishina asymp- 
totic form for isotropic photons in Eq. (|14|l is noticeable, high- 
lighting how in this domain, the extreme kinematic collimation 
imposed by the \ 2 > ^ 2 (1 — cos0)/2 condition is more con- 
straining than the beaming of the photon distribution within 



r]9 m > 1 



an angle 9 m . This characteristic is evinced in Fig. ^ via 
the approximate constancy of the ratio K,(a, r])/H(a, rj) for 
large 77 : comparison of Eq. (|23|l with Eq. (|14fl reveals that 
K(a, i])/H(a, 77) — + A(a)/(a + 1) as 77 — > 00 . 



p 
a" 

\ 
1? 
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Fig. 1. — The general behavior of K,(a, rj) is depicted as a function of rj 
as the heavy-weight curves, for three values of the spectral index a , as indi- 
cated. Also exhibited as the dotted cu rve is the asymptotic approximation 
for K(a, rj) for rj 3> 1 , precisely Eq. 1231 with the substitution 9 m — > 1 . 
To compare with the function Ti.(a, rj) that appears in the formalism for 
isotropic photons, in the upper portion of the plot the values of the ratio 
K,(a, 7j)/H.(a, rj) are illustrated as lighter weight curves, specifically for the 
indices a = 1, 2, 3 ranging from top to bottom. 

For the practitioner desiring an even more compact approxi- 
mation for the optical depth, a useful result that embodies the 
essence of both the two beamed forms in Eqs. (|17f) and 1231) 
can be put forward. There is no unique choice, however a com- 
parison of the coefficients in Eqs. I|16|) and (|18fl suggests the 
following mapping that mimics the appearance of the isotropic 
form: 



Q = min< 1 



X 



B(a) Wn 



(24) 



Insertion of this into Eq. Q then yields a development identical 
to that leading to Eqs. Hll|) and l|12l) . but with the substitution 
77 — > B(a) rj9 m . Hence, in cases where B(a) r)9 m < 1 , this 
protocol automatically generates the result in Eq. (|17fl . Com- 
parison of Eq. l(T4l . with 77 — * B(a)r)9 m , and Eq. then 
establishes the correspondence 



B(a) 



(a + l)A(a) 



l/(2a) 



(25) 



from these Klein-Nishina limiting forms, and the resulting com- 
pact approximation 



z. 



;(0m, V) « A{a)H(a, B(a) 77^) 9% 



(26) 



is generally accurate to around 10-20% for r]9 m J> 1 . Note 
that B(a) is of the order of magnitude of unity for typical a . 
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This approximation essentially maps over to a correspondence 
of B(a)8 m h 1/r between beaming angle and bulk Lorentz 
factor r of the emission region. Clearly such an approximation 
is subjective, but does not compromise the appropriateness of 
pair production attenuation calculations. Notwithstanding, for 
the investigation here, Eqs. (|21|l and (|22|1 will form the basis of 
the calculations below. 

2.3. General Constructs for Complicated GRB Spectra 

The adopted form for the approximate optical depth due to 
pair creation from a beamed conical power- law photon distribu- 
tion with a low energy cutoff at dimensionless energy s c (gen- 
erally <C 1 ) is 

r 77 (e) = n 7 a T R)c(a, 9 m ^feeZ) 9% e"" 1 . (27) 

Here, given that the angular form in Eq. Q is normalized to 
unity, the normalization parameter n 7 , defined in Eq. ( is 
such that n 1 9^/(a — 1) represents the total photon number 
density above the pair production threshold e ~ l/9 m for the 
beamed photon distribution. The expression corresponding to 
Eq. (|27|l for a power-law that extends down to zero energy but 
possesses a high-energy cutoff is obviously obtained by sub- 
tracting this result from that for an infinite power-law, namely 
Eq. (|T7|) . It is then a simple matter to construct optical depths 
for more complicated burst spectral structure. The particular 
spectrum that will be adopted throughout this paper is that 
used in Baring & Harding (1997a; hereafter BH97a), namely 
a power-law broken at a dimensionless energy £ B with a low 
energy cut-off at e c : 

rO, if e < e c , 

n(e) = n 7 e- ah < e^£" Qi , if £ c < e < £ B , (28) 
[e* h e- ah , £>£ B ■ 

This form mimics the generic shape of burst spectra (see Ta- 
ble 1 below for a listing of pertinent spectral parameters for 
BATSE/EGRET bursts), and the introduction of a sharp cutoff 
at £ c is to encompass scenarios where there is severe soft pho- 
ton paucity. Note that typical burst radiation mechanisms such 
as synchrotron emission and inverse Compton scattering from 
quasi-isotropic particles formally generate spectra with £ c — > , 
though in practice, the inverse Compton mechanism, so-called 
pitch angle synchrotron emission (e.g. see Epstein 1973), and 
moreover self-absorbed situations (e.g. Crider & Liang 1999; 
Granot, Piran & Sari 2000), produce extremely flat spectra so 
that finite choices of £ c can be realistic. Using the formalism 
described above, the optical depth for pair production attenua- 
tion of such a broken power-law/truncated photon distribution 
is found to be 

^g> « e{^,, c) -/cK, B )}^ 

(29) 

where 

rj B = 9 m ^/e B ~e , i] c = 9 m *Je^s . (30) 

Observe that at £ ~ 1/(e b 0^J the two spectral terms in Eq. (|2*9"|) 
become comparable, indicating that this energy is the domain 
where the attenuation "image" of the the BATSE-band spec- 
tral break at £ B is experienced. This expression for the optical 



depth, to be used hereafter in the computations of this paper, 
closely resembles that in Eq. (2) of BH97a, differing only by its 
use of K, functions in the place of TL functions and the replace- 
ment r — > 1 /9 m ; ensuing results in this paper are quantitatively 
but not qualitatively different from those in BH97a. 

Clearly, more gradual spectral curvature can be treated by fit- 
ting the GRB continuum with piecewise continuous power-laws, 
generalizing the structure inherent in Eq. I|29|l to sums of vari- 
ous power-law factors combined with differences in JC functions. 
The technique presented here is well suited to this adaptation, 
though for typical burst spectral forms, only modest changes 
in r 77 will be afforded by such refinement of the optical depth 
calculation, and then only deep into the attenuation trough. 
There may exist, of course, occasional pathological exceptions 
that require such refinement. 

2.4. Cosmological Background Targets 

The cosmic background infra-red (IR) and optical starlight 
fields provide target photons external to gamma-ray bursts that 
can attenuate in the EGRET/GLAST band. They effect linc- 
of-sight absorption, and so are exponential functions of the line- 
of-sight optical depth. Moreover, they are independent of burst 
characteristics such as flux, variability and spectral evolution, 
and therefore are markers only of the distance or redshift to 
a given burst. Such a property is touted as a means of prob- 
ing the background fields using any bursts (and active galaxies) 
detected by the atmospheric Cerenkov technique. The atten- 
uation effect of intervening cosmic background fields on hard 
gamma-rays via 77 — > e + e~ has been extensively studied in 
the context of active galaxies (e.g. Stecker, de Jager & Salamon 
1992, MacMinn & Primack 1996; for reviews, see Primack et 
al. 2001, Stecker 2001), with predictions being sensitive to the 
assumed level of the infra-red background, a historically poorly 
measured quantity. This uncertainty has been pervasive in the 
discussion of absorption of TeV photons from active galaxies, 
predominantly for AGNs that are detected at TeV energies by 
atmospheric Cerenkov telescopes have been nearby (e.g. Mrk 
421 and Mrk 501 are at z = 0.031 and z = 0.034 , respectively). 
However, this discussion has experienced a profound develop- 
ment in the recent observation (Aharonian et al. 2006) by the 
HESS telescope array of harder than expected TeV-band spec- 
tra in two more distant blazars, H 2356-309 at z = 0.165 and 
1ES 1101-232 at z = 0.186 (see also Albert et al. 2006b, for the 
very recent MAGIC detection of 1ES 1218+30.4 at 2 = 0.182 ). 
This extension to more distant sources has reduced the upper 
bounds on the extragalactic background IR light to interesting 
levels that are within 50% of lower limits set by galaxy counts 
with Hubble Space Telescope data. This advance moves the field 
of TeV gamma-ray astronomy much closer to being a powerful 
probe of these radiation backgrounds. 

Line-of-sight attenuation of GRB spectra provides a very dif- 
ferent case, at least for those bursts at moderate to high red- 
shifts that were seen by the BeppoSax mission, and are cur- 
rently being detected by Swift. Typical spectral attenuation 
expected in these sources was summarized by Mannheim, Hart- 
mann & Funk (1996), and impacts the sub-100 GeV band when 
z > 1 . Stecker & de Jager (1996) concluded that the 18 GeV 
photon seen by EGRET from GRB 940217 imposed a constraint 
of z <^ 2 for this source in order to evade interaction with in- 
tervening background photons. In the burst context, the higher 
redshift pushes the target photons into the near IR and optical 
starlight bands, which are better measured at low to moderate 
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redshift than the classical IR field. Additional uncertainty is en- 
tailed in the evolution of such backgrounds with redshift during 
the epoch of rapid star formation, yet galaxy count data can 
yield fairly well constrained models for the evolving background 
(Kneiske, Mannheim & Hartmann 2002). Since the focus here 
is on internal pair creation in bursts that is germane to the 
EGRET/GLAST band, calculations of attenuation by cosmo- 
logical backgrounds are beyond the scope of this presentation. 

3. ATTENUATION OF GRB SPECTRA: EVOLUTIONARY EFFECTS 

The GRB spectral forms adopted here have explicitly as- 
sumed, for simplicity, that there is only one radiation compo- 
nent contributing to the gamma-ray band. This is supported 
by the majority of EGRET burst data (see Table 1 for spec- 
tral indices), with the delayed 18 GeV photon in GRB 940217 
providing only a suggestion of a second component. Evidence 
for additional components is contingent upon broad spectral 
coverage. Accordingly an interesting possibility has been of- 
fered by the so-called MILAGRITO burst (Atkins et al. 2000; 
Atkins et al. 2003), GRB 970417, a BATSE burst that was 
seen with ~ 3<7 significance in the TeV band by the water tank 
Gerenkov detector MILAGRITO that served as the prototype 
for MILAGRO. The extrapolation of the BATSE spectral data 
up to the TeV band falls over four decades below the claimed 
MILAGRITO flux (see Fig. 9 of Atkins et al. 2003), thereby 
arguing strongly for the existence of a second component. No 
other bursts have been seen by the MILAGRITO /MILAGRO 
experiment, nor other TeV-band telescopes, with MILAGRO 
(Atkins et al. 2005) and MAGIC (Albert et al. 2006a) provid- 
ing interesting upper bounds to two bursts. 

Multi-component theoretical models exist in the literature 
(e.g. Meszaros & Rees 1994; Meszaros, Rees & Papathanas- 
siou 1994; Katz 1994; Dermer, Chiang & Mitman 2000), so 
the relationship between the current analysis and models with 
several radiation components cannot be summarily dismissed. 
Inspection of the developments of Section 2 quickly reveals that 
the pair production optical depth formalism presented here is 
entirely valid for multiple component spectra as long as they 
use the BATSE band component as the population of target 
photons. The super-GeV spectra may differ tremendously due 
to the presence of additional components, however the attenu- 
ation factors will depend only on the known BATSE/EGRET 
spectral forms. Accordingly, the signatures presented here will 
be widely applicable to burst spectra in the EGRET/GLAST 
band, unless additional components emerge below 100 MeV. 

In this paper, the analysis is restricted to cases of a single 
break in the single component power-law photon spectrum. This 
modest limitation is motivated by expediency, and can be relin- 
quished when more accurate representations of a GRB source 
spectrum are desired. Such improvements can be effected by 
approximating the spectrum by a piecewise continuous broken 
power-law with a sufficient number of components to achieve the 
desired accuracy. Yet such refinements rarely prove necessary if 
the specialized form adopted here in Eq. I|29|l is tailored to agree 
with published Band model fits (Band et al. 1993; Preece et al. 
2000). The reason for this centers on where the image of the 
target MeV photons is realized in the absorbed portion of the 
spectrum. The difference between the sharply broken power- 
law in Eq. I|28|l and the Band spectrum is manifested only near 
the break in the MeV band, leading to generally modest un- 
derpredictions of the attenuation factors, and then only near 
the depths of the troughs (that appear in Figures |2HH below) 
in the band above 10 GeV. This is the low-flux portion of the 



burst spectrum that will generally fall below current instrumen- 
tal sensitivities, and in some cases be obscured by attenuation of 
cosmological background radiation fields. Accordingly the sharp 
breaks assumed here shall suffice for the present investigation. 

3.1. Typical Absorption Characteristics 

To complete the determination of optical depths, one re- 
quires knowledge of the source emission region size R and 
the number density (or spectral normalization parameter) n 7 
of photons therein. For nearby bursts that sample Euclidean 
space, the coefficient n 7 can be obtained from the flux / 5 n 
at 511 keV per 511 keV energy interval via the simple re- 
lation /511 = Ann^cR 2 / d 2 , where d is the distance to the 
source. For such an absence of cosmological corrections, the 
perceived source size R can be written as R = R v r* , where 
r is the bulk Lorentz factor of the emission region with respect 
to the observer at infinity, and R v — cAt is the effective size 
of the region as deduced from observed variability timescales 
At. Following Baring & Harding (1996, 1997b), the value of 
R v = 3 x 10 7 cm is adopted here, commensurate with BATSE 
and COMPTEL variability timescales of At — ms, and an- 
ticipated variability timescales from the GLAST Gamma-Ray 
Burst Monitor (GBM). The index A is chosen according to the 
presumed geometry and structure of the emission region. If the 
constraining structure scale (e.g. of dumpiness, or between col- 
liding expanding shells) that defines source variability is along 
the line of sight (longitudinal) to the observer, then A = 2 . 
This is the preferred scenario in many GRB models that invoke 
colliding shells emanating perhaps from a central compact pow- 
erhouse. When the variability couples to structure transverse 
to the line of sight, then A = 1 , a situation that often appears 
in models of jet emission in active galactic nuclei. Such trans- 
verse variability was adopted by Baring & Harding (1997a, b). 
Theoretical arguments can be made for both A = 1 and A = 2 
scenarios (or hybrids thereof), though observational discrimina- 
tion between them is not yet afforded by the data. Both cases 
are considered in this paper, exhibiting only small qualitative 
differences between the two. 

The earlier analyses of Baring (1993, 1994) and Baring & 
Harding (1996, 1997a, b) mentioned cosmological corrections 
only in passing, and moreover in an incomplete fashion, pri- 
marily because these works predated BcppoSax burst detections 
that spawned redshift determinations via their afterglows. At 
this juncture, there is an observational mandate to include cor- 
rect redshift dependence in the pair production analysis. Such 
correction factors to the Euclidean (low z ) case enter in sev- 
eral ways. The simplest one connects to the source variability, 
whose intrinsic timescale is actually At/(l + z) . This is the only 
redshift contribution to the inferred source size if the Lorentz 
factor r is assumed independent of redshift. This restriction 
on r may be inaccurate, since there can be a significant (and 
as yet unknown) correlation between the nature of the burst 
progenitor and circumburst environment and its redshift. 

The major redshift dependence in the problem enters through 
the evolution of the source spectrum in its transmission to 
Earth. This receives contributions from the expansion of comov- 
ing volumes, which alters the perceived density of photons in the 
source, the redshifting of photon energy, and the frequency with 
which photons arrive at a detector. Spectral transmission cal- 
culations essentially connect to a luminosity transmission for- 
malism, which is textbook and can be found, for example, in 
Chapter 5 of Longair (1998). The core property is that the 
comoving, isotropic bolomctric source luminosity L bol and the 
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detected bolometric energy flux S hol are coupled via the lumi- 
nosity distance d L : 

O -^bol R Q foi \ 

bho1 - 4^ " 6grb ' z ■ (31) 

Here S'grb.z = £boi/(47ri? 2 ) is the GRB source flux at its outer 
periphery, i.e. at redshift z , and is a measure intrinsic to the 
source. While L bol is integrated over solid angles, Sboi is in- 
dependent of any collimation, provided the observer lies within 
the cone of collimation. 

The luminosity distance thus encapsulates the evolution of the 
energy fluxes in the curved spacetime of the evolving universe, 
essentially coupling to the conservation of energy in the comov- 
ing frame. Observe that both Sboi and S'grb.z possess units of 
energy per unit area per unit time (e.g. erg cm -2 sec -1 ). For a 
particular choice of cosmology, the form for d L is 

4 = WFV' (e) ' <32 » 

where Hq is Hubble's constant, and for k — sgn(f2 m + £1a — 1) , 

{sin O, k = 1, 

8, k = 0, (33) 

sinh 0, k = — 1, 

is the factor in the Robertson- Walker metric that describes the 
solid angle modification from flat space. Its argument is the 
development angle, and is given by 

Q = ^\T~9\J Z J^ , |l-fi|>0 , (34) 

with f2 = Q m + f^A being the total density expressed as a frac- 
tion of the closure density p c = 3Hq/(8itG) , and 

E{z) = Vn m (l + zf + fi A + (l - ft) (l + z) 2 . (35) 

Here, the radiation contribution to the universe's dynamics is 
assumed negligible, as it is for the epochs since burst genesis. 
This form for the luminosity distance is employed in standard 
expositions on supernova probes of cosmology (e.g. Perlmutter, 
et al. 1997; Riess, et al. 1998). 

In keeping with the results of those programs, and follow- 
ing on from the precise determinations of cosmology by WMAP 
microwave anisotropy observations (e.g. Spergel et al. 2003), 
here it is assumed that Hq = 0.72 km/sec Mpc -1 , and a choice 
of £l m — 0.27 and f^A = 0.73 is made (i.e. f2 — > 1), 
corresponding to a flat space, k = approximation. This 
simplifies the form of E{z) , and moreover the expression for 
<i L . Numerical integration of Eq. I|34|) is routine, and perti- 
nent asymptotic limits are w z — Cl\ for z <C 1 , and 

6 2^11- ^1/(1 - n A ) 1/3 F(l/6, 2/3; 7/6; 00 for z -> oo . 
These behaviors guide one to derive the useful analytic approx- 
imation for f^A = 1 — On = 0.73 : 

a ~ c z(l + z) \, z 2 \ 

dL ~ TT 1 + 0.29* \ 1_ 9(16T?) / ' ( j 

which is accurate to better than around 2% for < z < 40 , 
and generally better than 1% in the interval 1.7 < z < 40 . 

Eq. (|31|l can be easily converted into a spectral equivalent, 
and moreover one appropriate for photon number densities, by 



expressing the fluxes as integrals over differential photon spec- 
tra. For general spectral forms n(e) , the GRB source flux can 
be written as 

/>oo 

S'grb.z = 4:n m e c 3 J e z n(e z )de z 

= 4tt (1 + z) 2 m e c 3 J Sq n(jl + z] EqJ de , 

where e z = [1 + z]sq is the photon energy at the source cor- 
responding to an observed dimensionless energy of £q . An 
identical form to Eq. (|3*7jl is realized for 5 bo i , with the sub- 
stitution z — > , so that the relationship between the inte- 
grands of the two flux forms reproduces Eq. (7) of Woods & 
Loeb (1995). For the purposes of normalization in this paper, if 
an observer detects a differential photon number flux spectrum 
/( £ o) = /511 £0 Q above e c (with /511 in units of cm -2 sec -1 , 
and f5u/(a — 1) being the total number flux above 511 keV), 
assuming s c *C 1 for simplicity, then 

/>oo 

Shoi = fbiim e c 2 I 4 -a ffe , (38) 

J ec 

a result that is convergent for a > 2 . Then Eqs. I|31|) . 1371) 
and lEHl) can be combined to yield 

/ 5n = (0.511) 1-Q /(lMeV) = 4^n 7 c^ (1 + z) 2-a , (39) 

L 

independent of the choice of e c ■ Inversion of this result then 
leads to the specification of n 1 in terms of observables. Ac- 
cordingly, the redshift dependence of the optical depth r 77 is 
encapsulated entirely in the factor d 2 (l + z)" -2 , a result that 
can be inferred from Eq. (8) of Woods & Loeb (1995), and 
also Eq. (5.68) of Longair (1998). Note that this form con- 
trasts the erroneous claim made by Lithwick & Sari (2001) that 
r 77 oc (1 + z) 2 " -2 . In fact, just a single power a of 1 + z 
should appear since, due to pair creation threshold kinematics, 
the optical depth must trace one power of the inverse of the tar- 
get (i.e., observed) photon distribution, together with additional 
1 + z factors that are not directly connected to the spectrum. 

Note that since d L « cz /Hq = d for z <C 1 , Eq. i|3"9"|) reduces 
to the nearby source form /511 = 47m 7 cR 2 /d 2 in the z 1 
limit. In the implementations here, a — > ah is set. The rela- 
tionship between the theoretical flux parameter /511 and the 
the observed flux /(lMeV) at 1 MeV per MeV interval (tab- 
ulated for selected EGRET bursts in Table 1) is included in 
Eq. I|39|) for future reference. This concludes the formalism on 
cosmological input for the pair production attenuation problem. 

The attenuation of a theoretical gamma-ray burst continuum 
can now be computed. Neglecting the influence of pair cas- 
cading, if the emission region is comparatively confined, the 
attenuation of a spectrum like that in Eq. (|28|l is exponential 
in the optical depth r 77 given in Eq. (|29() . However, generally 
it can be distributed spatially with skin-depth effects in radia- 
tive transfer operating, so that the attenuation is by a factor of 
approximately 1/(1 + r 77 ) . Both cases are exhibited in the nu- 
merical results that follow. The principal spectral structure as- 
sociated with pair production attenuation is illustrated in Fig. [2 
for the BATSE/EGRET burst GRB930131, the so-called Super- 
bowl burst. This bright source exhibited the flattest spectrum of 
all EGRET bursts, and so is the prototype candidate for a burst 
easily detectable by GLAST at 100 MeV and above (see Band, 
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et al. 2004, Cohen- Tanugi, et al. 2004, or Omodei et al. 2006, 
for brief discussions of expected GLAST performance for burst 
detections). Since GRB 930131 pre-dated the determination of 
redshifts in optical afterglows, a redshift of z = 1 was assumed 
in the modeling (i.e., luminosity distance rf L = 6.54 Gpc), typi- 
cal of BeppoSax and Swift long-duration bursts. Relevant source 
data are listed in Table 1, which is a composite and extension of 
Table 1 of Baring & Harding (1997a) and Table 2 of Baring & 
Harding (1997b). While such flux-selected bursts might well be 
at redshifts lower than unity, GRB 990123 at redshift z > 1.61 
(Hjorth et al. 1999) provides a counterexample. 

TABLE 1 

Spectral Parameters for BATSE/EGRET Bursts 



BATSE" EGRET /COMPTEI/ 



GRB 


E B (MeV) 


oci 


Oth 


/(I MeV) 
(cm-'sec^MeV -1 ) 




(MeV) 


910503 


0.4 ± 0.2 


0.7 ±0.1 


2.1 


8.71 ± 0.49 


2.2 ±0.2 


170 


910601 


0.6 ±0.2 


1.0 ±0.1 


2.4' 


0.50 ±0.10'' 


2.8 ±0.2 i 


5.0 d 


910814 


1.2 ±0.1 


0.9 ±0.1 


2.3 s 


13.5 ±0.85 


2.8 ±0.2 


60 


930131 


0.7±0.1 


1.2 ±0.1 


2.5 


1.95 ±0.26 


2.0 ±0.2 


1000 


940217 


0.8 ±0.1 


1.2 ±0.2 


2.4' 


0.36 ± 0.03 


2.5 ±0.2 


3380' 


950425 


0.4 ±0.2' 


0.9 ±0.1' 


1.8 4 


1.62 ±0.09 


1.93 ±0.04 


120 


990123 


0.8 ±0.2 


0.6 ±0.1 


3.1 


1.7 ±0.2 


2.8±0.2 / 


18 



NOTE. — The time-averaged values of the break energy E B , the spectral index below ( ai ) 
and above ( ) the break, the (measured or extrapolated) source photon flux at 1 MeV 
per MeV bin (/(l MeV)), and the maximum energy e UAX of detection, for seven bursts 
observed with high significance by both BATSE and EGRET. All bursts listed were detected 
with the EGRET spark chamber, except for GRB 990123, which was seen in the TASC 
only (Briggs et al. 1999). (a) For the first three bursts, most BATSE parameters are taken 
from Band et al. (1993, the Band model fit), while BATSE data were obtained for GRB 
930131 from Kouveliotou et al. (1994), for GRB 940217 from Hurley et al. (1994), for GRB 
950425 from Catelli et al. (1996), and for GRB 990123 from Briggs et al. (1999). (b) These 
entries were obtained separately from an integration of the time-binned data archived in the 
electronic version of the BATSE Spectroscopy Catalog in Preece et al. (2000). (c) The 
high-energy spectral indices (o^), fluxes and maximum energies as determined by EGRET 
and/or COMPTEL are from Schneid et al. (1992: GRB 910503) and Kwok et al. (1993: 
GRBs 910601 and 910814), Sommer et al. (1994: GRB 930131), Hurley et al. (1994: GRB 
940217), Catelli etal. (1996: GRB 950425) and Briggs et al. (1999: GRB 990123). (d) These 
are COMPTEL measurements for GRB 910601 (see Hanlon, et al. 1994), quoted because 
COMPTEL detected it at a higher energy than did EGRET, (e) The notable 18 GeV photon 
from GRB 940217 was omitted here because it was not contemporaneous with any emission 
below 100 MeV. (f) This index from the OSSE/COMPTEL data was moderately consistent 
with the EGRET TASC value of 2.71 ±0.1 that was subject to poorer statistics (see Briggs 
et al. 1999). 

Spectra are presented for different Lorentz factors T = 1 /0 m 
(this coupling was adopted by Baring & Harding 1997a, though 
in general other T ~ l/8 m choices are equally acceptable), in- 
corporating only internal absorption and not treating propa- 
gational attenuation through background radiation fields. The 
spectra indicate marked absorption whose onset energy couples 
mostly to the photon flux, source distance and T , as detailed, 
for example, in the extensive investigations of Baring (1993) 
and Baring & Harding (1997b), and also depends weakly on the 
EGRET band spectral index ah - For 1/(1 + r 77 ) attenuation, 
above this turnover, the immediate spectral index is 2a/, — 1 , 
which flattens to ct^ + a; — 1 at energies higher than T 2 /e B , 
where the continuum below e B is sampled as the pool of tar- 
get photons. Such inferences can be made by inspection of the 
two terms contributing to Eq. I|29|) . These results resemble the 
character of the spectra illustrated in Baring & Harding (1997a) 
and Baring (2001), though note that higher Lorentz factors are 
appropriate here because of the greater luminosity distances in- 
volved (e.g. d L = 6.54 Gpc here as opposed to 1 Gpc in Baring 



2001). Spectra above around 30 GeV must, of course, be con- 
volved with external absorption due to the cosmic infra-red and 
optical backgrounds for typical burst redshifts. 
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Fig. 2. — The 7-7 attenuation, internal to the source, for GRB 930131 
for bulk Lorentz factors T = 500, 1000, 1500 of the emitting region, as- 
suming a redshift 2 = 1 typical of long duration bursts. The source 
spectrum (in e 2 /(e) format) was a power-law broken at = 0.7 MeV, 
with spectral indices a; = 1.2 and = 2.0 (see Table 1). Cases of 
1/(1 + t 77 ) attenuation (solid curves) and exponential attenuation (long 
dashed curve; T = 1000 only) are exhibited, for which a variability size 
of R = 3 X 10 7 r 2 cm was assumed. Also, the short dashed line depicts a 
r = 1500, 1/(1+t 77 ) attenuation example for _R = 3xl0 7 ITcm. The filled 
circle denotes the highest energy EGRET photon at ~ 1000 MeV (Som- 
mer et al. 1994). The threshold and sensitivity for the post-Tgo Whipple 
observations of later bursts (Connaughton et al. 1997) is indicated by the 
"WHIPPLE" box, while the MAGIC upper bounds for GRB 050713a (Al- 
bert et al. 2006a) that coincided in part with the prompt emission are 
denoted by the heavyweight histogram. 

The potential for observational diagnostics is immediately ap- 
parent. First, the extant EGRET data already provides a lower 
bound to r : the dot on Fig.[2]represents the highest energy pho- 
ton from GRB930131, and clearly suggests that T ^ 700 if the 
burst was at a redshift of z ;> 1 . It is also clear from the figure 
that the exponential and 1/(1 + r 77 ) attenuation signatures are 
strikingly different, so that a sensitive hard gamma-ray detector 
will easily be able to discriminate between these two possibili- 
ties. While realization of the 1/(1 + r 77 ) case is anticipated, 
should exponential attenuation be observed, it would be a pro- 
found diagnostic indicating that the emission region is relatively 
small, confined to a scale length shorter than the mean free path 
for pair attenuation. Comparison of the R = 3 x 10 7 L 2 cm and 
R = 3 x 10 7 F cm examples for F = 1500 in Fig. |3 indicates 
that, in the absence of time- varying spectra, it will be difficult 
to discriminate between these longitudinal and transverse vari- 
ability cases, since they cannot be clearly resolved from Lorentz 
factor inferences. Specifically, R cx T scenarios require slightly 
higher T for transparency at a given photon energy than do 
longitudinal variability cases. 

To aid in the interpretation of pair attenuation it is useful 
to analytically determine the mean interacting (target) photon 
energy s iat (s) for a given observed photon energy e. This is 
computed as the ratio of an interaction probability like that 
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in Eq. (0), but with the extra weighting factor u> in the inte- 
grand, to the interaction probability with unit weighting factor. 
Retracing the analytic development in Sec. 12.11 it is quickly de- 
duced, using the spectral form in Eq. (0} for the target photon 
population, that 



£int(£) 



4 /C(a - 1, rjc9m) 
/C(a, r] c 9 m ) 



(40) 



For this spectral specialization, e c plays the approximate role 
of the break energy e B , though an alternative and slightly more 
involved form can easily be obtained to pertain to the spectrum 
in Eq. (|28l) . Ratios of JC functions can be inferred from Fig.^ 
and are typically of the order of a few for r\ ~ 1 . Accordingly, 
for the specific example in Fig. [21 Eq. (|40() clearly establishes 
that the attenuation image of the BATSE-band spectral break 
at e B appears at the energy e ~ 4r 2 /e B , which corresponds 
to the pair threshold condition in the bulk fluid rest frame, i.e., 
(e/r) (e B /r) ^> 4. This is generally at or above the threshold 

energy of atmospheric Cerenkov telescopes. 

The first intensive TeV-band exploration of bursts was the 
Whipple Telescope campaign (Connaughton et al. 1997), whose 
results are typified by the "WHIPPLE" box in Fig-HJ This cam- 
paign in 1994-1995 consisted of targeted observations of nine 
BATSE bursts within, at best, around two minutes of trigger, 
coordinated via the BACODINE network. For each burst, the 
Whipple observations comprised 6 scans of around 30 minutes 
each with slightly offset sky positions (by w 3° ) to cover the 
BATSE error boxes. No TeV-band detections were made, and 
so the WHIPPLE box in Fig. [2] represents upper limits to the 
potential flux from the bursts studied in Connaughton et al. 
(1997). Due to the slewing times of a few minutes, these bounds 
are more germane to the burst afterglow phase, for which the 
MeV-band flux level is unkown. In contrast, the recent rapid- 
slewing observation (Albert et al. 2006a) of GRB 050713a by the 
MAGIC Telescope, triggered by a Swift alert, enabled a viewing 
of the burst at about 40 seconds after its onset, for a period of 37 
minutes (Swift determined Tgo = 70 ± 10 sec). No gamma-ray 
signal was detected by MAGIC, and the resulting upper bounds 
are indicated in Fig. |21 as the heavyweight histogram. Since 
no redshift was determined for GRB 050713a, this sensitivity 
histogram corresponds to the z = 1 flux upper limit listed in 
Table 1 of Albert et al. (2006a). This flux limit is approxi- 
mately commensurate with the fluence limits obtained by the 
MILAGRO Cerenkov water tank experiment for GRB 010921 
(Atkins et al. 2005), the only burst surveyed in the MILAGRO 
campaign of 2000-2001 that possessed a moderately low, mea- 
sured redshift ( z = 0.45 ). STACEE also obtained a flux upper 
limit to GRB 050607 (Jarvis et al. 2005). The MAGIC sensitiv- 
ity depicted in Fig. [3 is representative of the present capability 
of state-of-the-art atmospheric Cerenkov systems in terms of 
prompt GRB observations. 

3.2. Time Evolution of Spectra 

The photon counting statistics associated with EGRET burst 
detections generally only permitted using time-integrated spec- 
tra for the purposes of pair attenuation considerations. The en- 
hanced sensitivity offered by GLAST by late 2007 may permit 
some exploration of time-dependent attenuation characteristics 
that are previewed in this Section, particularly if they emerge 
in bright bursts in the energy band below around 1 GeV. 

There are a plethora of time-dependent gamma-ray burst 
spectral models in the literature, often pertaining to afterglow 



evolution. For the prompt phase, there is evidence of preva- 
lent hard-to-soft evolution on average in the BATSE database, 
though the spectral hardness development during a burst can 
often be more complicated (e.g. Ford et al. 1995, Preece et al. 
1995). If the Lorentz factor T is effectively constant, then this 
observed variation could be due to evolving dissipational char- 
acteristics in the internal shocks, which are subtle and difficult 
to encapsulate in a cohesive manner. As an alternative scenario, 
the spectral evolution could primarily reflect the change in T on 
a timescale of the burst duration. For long-duration bursts with 
Tgo ~ 1-100 sec, the physical scale associated with their dura- 
tion is r 2 cT 90 ~ 3xl0 15 -3xl0 18 cmfor T ~ 300-1000. Hence 
bulk deceleration during the burst prompt phase is likely, and is 
a natural driver for GRB dissipation (e.g., Rees and Meszaros 
1992). Both evolutionary possibilities are examined here. 

The focus is first on a constant T , constant ah scenario with 
the flux allowed to vary, for which representative spectra are 
displayed in Fig. H for the case of GRB 950425. There, the 
emission spectrum evolves only by an across-the-board decline 
of flux with time, displayed in increments of half a decade per 
"snapshot" (solid curves). Accordingly, this case does not ex- 
plore adiabatic deceleration issues, or the role of radiative cool- 
ing, but does provide a definite correlation between evolution 
of the BATSE-band flux and the turnover energy E t in the 
GLAST window. To elucidate this interesting diagnostic, ob- 
serve that since the Lorentz factor is held constant, and at the 
turnover energy the last term in Eq. (|29|l provides the dominant 



contribution to t 77 , one finds that r 77 



cx f(l MeV)E? h - 1 



From this, it follows that the turnover energy couples to the 
MeV-band flux according to 



E t oc [/(l MeV) 



|-V(«fc-i) 



(41) 



established using a r 77 ~ 1 criterion. For the GRB 950425 ex- 
ample in the Figure, this gives E t ~ [f(l MeV)] -1 / 2 . Com- 
bining this result with the /(l MeV) E 2 ~ ah spectral form far 
above the break energy Eb then quickly yields 

E 2 t f{E t ) cx E?- 2ah cx [/(lMeV)]^^^- 11 (42) 

as the mathematical form for the locus of turnover points in 



3 — 2 a?, 



E, 



-0.92 



vF v space. For GRB 950425 this generates E\ 
since cth — 1.93, which is represented by the dashed curve in 
Fig. El Clearly, this correlation is extremely pronounced, and 
is a very useful probe of the constancy of T if the high energy 
spectral index does not evolve. In bright, flat-spectrum sources 
such as GRB 950425, GLAST may well be suited to investi- 
gating whether such character is exhibited, provided that the 
turnovers appear at below 1 GeV so that the LAT can accumu- 
late enough photons in a typical burst duration. In particular, 
for the r adopted in this example, clearly the relative probabil- 
ity of seeing a super-100 MeV photon to observing a 1 MeV pho- 
ton is enhanced at later times; whether this is connected to the 
detection (see Hurley et al. 1994) of the 18 GeV photon in GRB 
940217 over an hour after burst trigger is undetermined. Note 
also that the attenuation results apply to any time-ordering of 
the curves, not just monotonically declining fluxes. Accordingly, 
one would infer a decline of the turnover energy in time in a ris- 
ing flux precursor to the prompt GRB emission. 

Another property of the attenuation is that well above the 
turnover energy, the results are approximately independent of 
/(l MeV) , when e B is held fixed in time. This arises in the 
"recovery domain," where the target photons are below e B and 
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the optical depth r 77 is declining with energy till eventually the 
internal attenuation becomes insignificant (illustrated in Baring 
& Harding 1997b). At such energies higher than r 2 /e B , as 
oti = 0.9 for GRB 950425, the optical depth is almost indepen- 
dent of the photon energy so that it mostly scales as /(l MeV) . 
This behavior almost exactly cancels the normalization of the 
emission continuum so that the emergent spectrum is virtually 
independent of the MeV-band flux. For this source, given the 
location of the ACT sensitivity indicator in Fig. [3] (i.e. MAGIC 
flux upper limits for GRB 050713a), and the probable role of 
additional line-of-sight attenuation due to cosmic background 
radiation fields, the recovery domain is unlikely to be accessible 
to current generation Cerenkov telescopes. 
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Fig. 3. — The 7-7 attenuation, internal to the source, for GRB 950425 
for bulk Lorcntz factor V = 500 , assuming a redshift z = 1 typical of long 
duration bursts. The source spectrum (in e 2 f(e) format) was a power-law 
broken at i?B = 0.4 MeV (solid curves), with spectral indices Q; = 0.9 and 
ah = 1.93 (see Table 1). Only 1/(1 + T77) attenuation is depicted, and a 
variability size of R = 3 X 10 7 T 2 cm was assumed. The peak flux evolution 
between the curves is by successive factors of vTo from top to bottom. 
The dotted curve illustrates a hard-to-soft evolution situation and is vir- 
tually coincident with the corresponding fixed Eb = 0.4 MeV case above 
the break, and below 500 GeV. The filled circle denotes the highest energy 
EGRET photon at ~ 120 MeV (see Table 1). The attenuation turnover 
energy Et evolves, tracing out an E t -0 ' 92 locus, as discussed in the text. 
As in Fig.|2] the MAGIC flux limit upper bounds for GRB 050713a (Albert 
ct al. 2006a) are denoted by the heavyweight histogram. 

Finally, for this example, observe that since e B is not influ- 
ential in determining the turnover energy, it was held constant 
for illustrative purposes. In Fig. |3| the dotted curve depicts a 
case where strong hard-to-soft evolution would arise subsequent 
to the peak flux snapshot. Well above e B , and below the re- 
covery domain at E > 1 TeV, this curve is not distinguishable 
from the corresponding solid curve for the fixed e B case, since 
the target energies for photons near the turnover are above the 
break energy. In the E ;> 1 TeV portion of the spectrum, the 
extra supply of target photons below 0.4 MeV relative to the 
e B = 0.4 MeV solid line case with the same /(l MeV) clearly 
enhances the pair opacity and lowers the spectrum accordingly. 
Fig-Elencapsulates high energy spectral attenuation characteris- 
tics appropriate for constant T scenarios that admit evolution of 



the break energy e B with time, and therefore possesses broader 
applicability. The import of this broad scope will soon become 
apparent, when scenarios incorporating T evolution with time 
are addressed just below. Observe also that while not illus- 
trated in this example, evolution of ah with time can be easily 
accommodated to provide useful diagnostics on T . 

The simplest (but by no means uniquely preferred) case cor- 
responding to appreciable reductions of T on a timescale of 
the burst duration is defined by the hydrodynamic sweep-up by 
a blast wave of exterior material from the interstellar medium 
or a progenitor wind. This is essentially the scenario explored 
by Rees and Meszaros (1992), and encapsulated nicely in the 
paper by Dermer, Chiang & Bottcher (1999). It is this lat- 
ter analysis that is employed here, specifically the self-similar 
solution for blast wave deceleration in a uniform or spatially- 
structured circumburst environment. Assuming, for simplicity, 
that the density of material the blast wave impacts is uniform, 
then the sweep-up is precisely a relativistic extension of the 
classic Sedov- Taylor phase solution for supernova remnant ex- 
pansion. This hydrodynamic problem for an expansion with 
an ultra-relativistic equation of state was treated by Blandford 
& McKee (1976). In the limit that the magnetic field plays a 
passive role, the conservation of energy and momentum in an 
adiabatically-expanding, non-radiative scenario yields a radial 
form for the deceleration of 

r(r) cx r- 3 / 2 oc i- 3 / 8 , t dyn < t < t rad (43) 

after the coasting phase has passed on a timescale of t dyn . If the 
expansion is radiative, i.e. loses energy via radiation on times 
£ rad shorter than dynamical timescales t dyn , then the evolution 
is described by a momentum-conserving "snowplow" form (e.g. 
Blandford & McKee 1976; Dermer, Chiang & Bottcher 1999) 

T(r) oc r~ 3 oc i~ 3/7 , i rad < t dyn < t . (44) 

The dynamical sweep-up scale ct dyn is directly expressible in 
terms of the "explosion" energy £ , the solid angle Slo of col- 
limated expansion, coasting phase (initial) Lorentz factor Tq , 
and mass density p of the circumburst medium as 

*- - ts {ojt?r ■ (45 » 

Here the presence of the Tq factor out the front accounts 
for relativistic time dilation between the blast and observer's 
frames (similar factors mediate the r-t coupling in Eqs. |43j 
and 02); an d the other Tq factor couples to the energy conser- 
vation during sweep-up; these factors provide the only relativis- 
tic modification from the classic Sedov- Taylor phase timescale. 
Other radial/temporal dependences for T can be envisaged (e.g. 
Meszaros, Rees & Wijers 1998), for example that incurred by a 
power-law radial dependence for the circumburst density profile. 

The forms in Eqs. I)43|) and (|44|) suffice for the purposes of this 
exposition, yielding time-dependence templates for scoping out 
pair creation attenuation signatures. Collimation of the outflow 
causally impacts the kinematic phase space only very late in the 
GRB afterglow. The other ingredients needed are the Lorentz 
factor dependences for the spectral break energy e B and the flux 
at this or some fixed higher energy. These quantities are model- 
dependent. Adhering to the prevailing paradigm that the GRB 
emission mechanism is synchrotron radiation, Dermer, Chiang 
& Bottcher (1999) argued in their external shock model that 
the break energy, if it corresponds to synchrotron photons from 
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electrons of the minimum Lorentz factor "f n 
distribution, would behave according to 

£b oc [r(r)] 4 . 



in a non-thermal 



(46) 



The origin of this form assumes that 7 min oc T arises during en- 
ergy transfer in the sweep-up phase, that the ambient magnetic 
field is boosted by T in the comoving frame, and also incorpo- 
rates one power of T for the energy blueshift to the observer's 
frame. All these are reasonable contentions. Yet, it is clear that 
a sizeable fraction of bursts cannot generate low energy spec- 
tra consistent with the synchrotron mechanism using a cutoff 
distribution; this is the so-called "line of death" issue raised by 
Preece et al. (1998, 2000), discussed also in Baring & Braby 
(2004). Hence, the possibility of £ B (r) dependences other than 
in Eq. (|46|l can also be entertained. 

The flux evolution with time or Lorentz factor is encapsulated 
in the third part of Eq. (20) of Dermer, Chiang & Bottcher 
(1999). Physically, this form assumes that there is prompt dis- 
sipation of a sizeable fraction the available energy, namely the 
energy of matter swept up in the deceleration phase. For en- 
vironments where the the circumburst medium is uniform, this 
form yields a temporal dependence for e 3> £ B of 



where 



l(e) OC t-3v/(8-r )r ) K £ -a h TV 



A 16 , 2 

<p = Aa h - -g- + g rj r 



(47) 



(48) 



is the index defining the bulk Lorentz factor dependence for the 
source flux, and 



Vr = 



, adiabatic, 

1 , radiative 



(49) 



is the parameter that demarcates radiative and adiabatic expan- 
sion cases. Observe that with this definition, T oc t~ 3 /( 8 ~')>-) 
describes the cases in Eqs. ()43[) and ()44|) . so that the identity 
g = 3/(2 — rj r ) provides a connection to the parameter g that 
Dermer, Chiang & Bottcher (1999) use throughout their inves- 
tigation, with r oc r~ 9 . Observe also that inserting Eq. (14611 
into Eq. (gTJl yields 



ei n 



(e B ) oc t -i-3W(8-»h.) oc r 2 ( 4 +"")/ 3 



(50) 



commensurate with the evolution of the peak flux that is given 
in Eq. (16) of Dermer, Chiang & Bottcher (1999). Adiabatic 
expansions generate e| n(e B ) oc oc r 8 / 3 , a form dictated by 
energy /momentum conservation in sweep- up, which because of 
scale-invariance is deducible directly from Eq. 145|) ; radiative ex- 
pansions produce e^n(e B ) oc £~ 10/7 oc r 10 / 3 . The flux evolution 
prescribed by Eq. ljT7l) can be easily modified to accommodate 
expansions into non-uniform circumburst media, yet suffices for 
the illustrative purposes of this paper. 

It should be noted that while the temporal evolution employed 
in Dermer, Chiang & Bottcher (1999) envisages an external 
shock model for the prompt emission, and can directly be ap- 
plied to afterglow considerations, it is also germane for internal 
shock models if the gamma-ray activity arises during a deceler- 
ation epoch. Then, the colliding shocks would possess different 
Lorentz factors, which would on average slow down with the ex- 
pansion. The relevant T for the purposes of this analysis would 
then be some average of those of the colliding shells, and would 
be controlled by the hydrodynamics discussed just above. 
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Fig. 4. — An evolving sequence depicting 7-7 attenuation, internal to the 
source, for GRB 910503, again assuming a redshift 2 = 1 typical of long 
duration bursts. The source spectrum (in e 2 /(e) format) was a power-law 
broken initially at Eb = 0.4 MeV (uppermost solid and dashed curves that 
are coincident), with spectral indices aj = 0.7 and Q/, = 2.2 (see Table 1) 
that are maintained throughout. For this initial case, the bulk Lorentz 
factor was T = 600 . The highest energy EGRET photon at ~ 170 MeV 
(Schneid et al. 1992) is denoted by the filled circle. Again only 1/(1 + t 77 ) 
attenuation is depicted, and a variability size of R = 3 X 10 7 r 2 cm was 
assumed. The peak flux and hard-to-soft evolution between the curves is 
indicated via declining Lorentz factors, i.e., cases T = 365 and T = 220 as 
labelled. Th e hard-to-soft evolution of the break en ergy eb is prescribed 
via Eq. 1461 . while the flux evolution is given in Eq. 1471 . and is discussed 
in the text. Solid curves correspond to adiabatic blast wave deceleration, 
for whic h the attenuation turnover energy approximately traces out [see 
Eq. 1521 1 the depicted E\ 2 locus, while dashed curves illustrate a strongly 
radiative case, for which the locus of -r-v 



tratcd E 



2.0 



1 turnover points is the illus- 



line. 



Representative spectra for evolving F scenarios are exhib- 
ited in Fig. H for the case of GRB 910503. This elaborates 
upon the evolutionary sequence addressed in Baring (2001) for 
GRB 930131. Here ai and ah are held constant, so as to sim- 
ply demonstrate the key characteristics, though for some bursts 
(e.g. see Gonzalez et al. 2004) there are significant variations 
of ah with time. The flux and break energy evolve according 
to Eqs. H46f> and 1471) respectively, and in this example e c = . 
Notably, the turnover energies are now declining functions of 
time, contrasting the constant T example in Fig. |3 for which 
the turnover locus is E^~ 2ah . This behavior is obviously caused 
by the increased opacity due to the decline of T , which more 
than offsets the tendency of a reduced flux at typical super- 
MeV target photon energies to reduce the pair opacity. Such a 
signature renders it easy to discriminate GRB environs where 
bulk deceleration is prevalent during the prompt phase, from 
the case of constant T , where internal shock dissipation is com- 
pletely decoupled from external shock evolutionary dynamics. 
This is a powerful observational diagnostic that is only modi- 
fied in detail, but not in character, by model nuances such as 
spatial non-uniformity in the circumburst medium. 

A comparison of the solid and dashed curves in Fig. 0] indi- 
cates that it will be harder to discern between adiabatic and 
radiative evolutionary scenarios, though it may be marginally 
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possible. The tracks of the spectrum around the break energy 
do not distinguish substantially between the rj r = and 77,. = 1 
cases in that the separation of fluxes for given break energy is 
only modest, and break smoothing could render such trends in- 
distinguishable. The turnover feature adds extra information 
that is helpful, provided that the hard gamma-ray continuum 
is detected with sufficient significance over around at least 2 
decades of flux evolution at above 1 GeV. The Figure marks 
two loci that the turnovers trace in time for the respective adi- 
abatic and radiative cases. As for the constant T case, these 
loci define the correlation between evolution of the BATSE-band 
flux and the turnover energy E t , which can again be found us- 
ing the last term in Eq. (|29|) as the dominant contribution to 
r 77 . One finds that r 77 cx /(l MeV) E"^ 1 /T x+2ah , where 
A = 1,2 is the index defining the Lorentz contraction relating 
the perceived source size to the variability timescale. From this, 
and Eq. I|47|) . it follows that the turnover energy satisfies 



Et 



^X+2a h -<p)/(a h -l) 



(51) 



expressing the t 77 ~ 1 criterion, with ip as defined in Eq. I|48l) . 
Inverting this result and combining it with Eq. 147|) . the time- 
dependent spectral form valid far above the break energy e B , 
results in 



E 2 f(E t ) cx El 



g = 2 - a h 



tp(a h - 1) 



A + 2ai 



(52) 



as the mathematical form for the locus of turnover points in 
vF v space. For GRB 910503 with ah — 2.2 , and assuming 
A = 2 , this approximately generates E\ 2 for the adiabatic case 
(rj r =0) and Ef ° for radiative regimes (r/ r — 1), which are 
represented by the lightweight solid and dashed lines in Fig. 0] 
These are clearly distinct, and offer the attractive possibility 
of discriminating between adiabatic and radiative evolutionary 
scenarios for a given burst, particularly if explored in conjunc- 
tion with the e B evolution. 

Since g is a rapidly increasing function of ah , flat-spectrum 
bursts with high super-100 MeV fluxes offer the most promising 
discrimination: ah = 2 yields g — 4/5 for adiabatic (r] r = Q) 
regimes, and g = 5/4 for radiative ones (r] r = 1). However, 
note that when adopting a variability size of R = 3 x 10 7 r cm 
(i.e., A = 1), the loci in Fig. 21 change significantly, to ap- 
proximately Ef and Ef- 7 for adiabatic and radiative cases, 
respectively. Hence, in practice it may prove difficult to cleanly 
distinguish between the effects of non-adiabaticity in the ex- 
pansion and the nature of structure in the emission region using 
the turnover energy as a diagnostic. Using the sub-MeV band 
temporal traces of the break energy in conjunction should help 
separate these competing influences. 

4. THE POTENTIAL FOR OBSERVATIONAL DIAGNOSTICS 

As a benchmark for the GRB phase space that is germane 
to internal 77 opacity diagnostics by GLAST and perhaps at- 
mospheric Cerenkov telescopes, and future experiments down 
the line, one can map the contours the turnover energies trace 
for burst parameters typical of the bright sources EGRET that 
detected. For given T and ah , the maximum photon energy 
E MAyL — E t (in MeV) attainable in the EGRET-band power- 
law, before pair attenuation turnovers appear, is simply deter- 
mined by noting that the spectral structure at and below the 
BATSE-band break is immaterial to its evaluation. Therefore, 
one can set e B — > in the analysis, and then impose a r 77 ~ 1 



criterion to establish, using Eqs. © and l(T7|l . that 
511 T 2ah ^ ah ~ 1 ^ 



E* 



{n 7 cr T i? A(a h )H(a hl 1)| 



MeV 



(53) 



This energy is plotted as a function of for different T and 
fluxes /(l MeV) in Fig. El specifically for source luminosity dis- 
tances of 6.54 Gpc (i.e. 2 = 1). Such curves were determined 
assuming a variability size of R = 3 x 10 7 r 2 cm , and are appli- 
cable for both 1/(1 + t 77 ) attenuation and exponential attenua- 
tion. This depiction represents an update of Figure 6 of Baring 
and Harding (1997b) that reflects the moderate redshifts now 
associated with GRBs in the Swift/HETE era; Jakobsson et al. 
(2006) report a mean redshift of z ~ 2.8 in their compilation 
of 16 Swift bursts with measured redshifts. 
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Fig. 5. — The phase space diagram for the observational parameters 
Emax and , consisting of contours of c onstant T , as labelled. The 
contours correspond to solutions in Eq. 1531 of T- n = 1 . Here z = 1 , and 
a variability size of ij = 3 X 10 7 r 2 cm was adopted. The solid and dashed 
( r = 500 case only) contours correspond to the indicated choices of the 
flux /(l MeV) at 1 MeV. The observed EGRET values of E M ax and a h 
for the seven bursts in Table 1 are plotted as points; these bursts all have 
values of /(l MeV) somewhat different from the chosen contour values. 

The £ MAX curves are monotonically increasing functions of 
ah and T , and are decreasing functions of /(l MeV) , as indi- 
cated by comparison of the solid and dashed T = 500 examples. 
Since the optical depth is fixed at unity at the turnover energy, 
and the flux is pinned at 1 MeV, thereby providing a "pivot 
point" in the spectrum, increasing ah then lowers the number 
of interacting photons above pair threshold (generally above 1 
MeV in the observer's frame) so that E MAX correspondingly in- 
creases to compensate. The depiction of (ah, E MAX ) positions 
for the EGRET bursts listed in Table 1 clearly indicates that 
r > 100 in their emission regions if they are all at or near a red- 
shift of z = 1 . Observations of further extension of their power- 
laws by GLAST could push these Lorentz factor bounds consid- 
erably higher: much of the phase space below E MAX ~ 10 4 MeV 
and ah < 2.5 is readily accessible to GLAST, underlining its 
potential as a tool for probing GRB bulk motion characteristics. 
As the hard gamma-ray spectral index increases to ah ^> 2.8 , 
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the total photon counts collected by the LAT instrument will 
correspondingly be reduced, lowering the maximum energies to 
which it can probe to below 1 GeV. 

The results of Section 13.21 highlight one and possibly two 
groundbreaking probes of the burst environment by hard 
gamma-ray experiments that explore the GLAST window. 
First, the comparison of Figs. |3] and 0] offers strikingly differ- 
ent evolution of the turnover energy E t with time (i.e., declin- 
ing MeV-band flux) for constant T and evolving T scenarios. 
A sensitive hard gamma-ray detector would easily be able to 
discern between these two signatures, using data binned into 
several time windows, and cast light on the T evolution for a 
number of bright, energetic bursts. Moreover, these internal 
absorption signatures can be easily distinguished from external 
absorption characteristics due to background fields, which gen- 
erate turnover energies independent of time or MeV-band flux. 

This tracking of T evolution is a global possibility that is not 
confused by specific model nuances such as circumburst non- 
uniformity and departures from temporal monotonicity for the 
flux. Should a burst exhibit T = const, behavior, it can be 
concluded that the expansion timescale is shorter than the dy- 
namical one, so that the burst duration T90 in conjunction with 
Eq. I|45|) can be used to provide an approximate lower bound to 
the explosion energy £0/^0 per unit solid angle of the expan- 
sion, modulo the density of the environment surrounding the 
burst: 

^ ;> U = irgpc 5 T 9 3 , r = r = const. (54) 

On the other hand, should T changes be observed, whether 
increasing or declining with time, one could argue that the dy- 
namical timescale is smaller than the expansion age and that 
r < To . Accordingly, an upper bound £0/^0 < Uq on the 
collimated explosion energy would be obtained, and this could 
be tightened to an approximate estimate £0/^0 ~ T s pc 5 Tg , 
with considerable uncertainty due to the spread in T estimates 
derived from the pair attenuation turnovers. This calibration 
of £0/^0 is analogous to inferences made in non-relativistic 
supernova remnant expansions either in the coasting or Sedov 
phase. Furthermore, if a radiative regime is deduced to be oper- 
ating, then hard gamma-ray observations can also constrain per- 
tinent environmental parameters for a chosen radiation mecha- 
nism, via bounds on the cooling timescale. Notwithstanding, it 
is noted that marked variations of with time can degrade 
this diagnostic possibility. Such a property is evinced in spec- 
troscopy analyses of some bursts (see, for example, the joint 
BATSE/EGRET TASC study of Gonzalez et al. 2004), though 
the average characteristics of time-resolved ah determinations 
by BATSE are similar to the time-integrated ah values (e.g. 
Kaneko et al. 2006). GLAST will solidify the understanding of 
temporal variations, or otherwise, in ah for many bursts. 

An interesting, though somewhat less promising diagnos- 
tic is the discrimination between adiabatic and radiative/non- 
adiabatic evolutionary scenarios. While the evolutionary loci of 
E t for these cases are distinct in Fig. they are only modestly 
so, and can easily be blurred by breadth of the turnover, the 
influence of non-uniformities in the burst environs, or compli- 
cated variability in the spectral shape with time. Furthermore, 
significant changes in the value of the index g in Eq. (|52|l are 
incurred in going from longitudinal ( A = 2 ) to transverse vari- 
ability ( A = 1 ) situations. Hence it becomes harder to disentan- 
gle a variety of effects given an observed, evolving, GeV-band 
spectrum. Degeneracies can be partially removed by exploring 



correlations with the MeV-band peak evolution, in both flux 
and break energy. For example, Eqs. H46[) and H5U|) combine to 
give a locus of breaks whose slope depends only on rj r , and is 
independent of ah ■ How effective such a correlated diagnostic 
will be will depend on the sharpness of the respective breaks. 

GLAST's Large Area Telescope (LAT) offers potential for 
exploring these signatures, being considerably more sensitive 
than EGRET; its principal capability will be to provide lower 
bounds to T via its refinement of E MAX , and perhaps also 
the first observations of turnovers in the brightest subset of 
the burst population (e.g. see Omodei et al. 2006, for simu- 
lated spectra with detectable turnovers). Band, et al. (2004) 
and Cohen- Tanugi, et al. (2004) provide brief outlines of 
the anticipated performance of GLAST for burst detections. 
Being ephemeral in nature, GRBs are vastly different from 
steady sources, for which the LAT Performance Web Page 
http : //www-glast . slac . Stanford . edu/ sof tware/IS/glast 
_lat_performance.htm offers putative integral sensitivities for 
the 5cr detection of a high latitude source above a given energy, 
resulting from a one-year, all-sky survey. With the average 2 
steradian field of view for the LAT, this corresponds to around 
t lyl = 5 x 10 6 seconds observing a given source, far above the 
typical duration of bursts. For steady sources, background is 
an issue (especially in the galactic plane: see Hunter et al. 
1997, for the presentation of EGRET background determina- 
tions), but its influence on the observability for bursts virtually 
disappears since they are of such short duration. Gamma-ray 
burst detectability by the LAT correlates to the total number 
of photons counted above its threshold energy of 30 MeV. Po- 
tential instrumental deadtime due to data readout after triggers 
is small, about 5%. LAT Team simulations of detector sensi- 
tivities (Omodei, et. al. 2006) indicate that only a handful of 
bursts seen by GLAST each year will provide several hundred 
photon counts in the LAT above the threshold of 30 MeV, and 
perhaps one a year with of the order of a 1000 or so (and around 
one a year with 100 photons above 1 GeV). Such cases could, 
in principal, yield of the order of 30-100 photons above 300 
MeV if the burst possesses a flat E~ 2 spectrum. Observing 
turnovers around 100 MeV-T GeV would then be viable with 
such a database, and moreover, it may prove possible to track 
spectral evolution at around 100 MeV-1 GeV by binning data in 
two or three time intervals, with corresponding flux reductions 
by a factor of 2-5 in the MeV band (or at higher energies). 

Such is the optimal situation. Nature, of course, may not pro- 
vide turnovers at such low energies, so GLAST's potential for 
casting light on the above bulk Lorentz factor diagnostics will 
only truly be known after launch. Yet, bright EGRET bursts 
typically have flat spectra, on average around E~ 2 (see Din- 
gus, 1995, for a compilation encompassing most of the bursts 
considered in this paper), which may or may not be the result 
of flux selection effects at higher energies. These are ideal can- 
didates for exploring spectral turnovers in the super-GeV band, 
which were simulated by the LAT Collaboration in Omodei et al. 
(2006). At best, only an elite handful of bursts over the GLAST 
lifetime might be amenable to observing both turnovers in the 
sub-GeV band, and partially tracking their evolution with time. 
This number will increase if the bulk Lorentz factors are ac- 
tually slightly lower than those used in the illustrations here. 
Accordingly, GLAST might be capable of performing some of 
the GRB source diagnostics discussed in the text. In addition, 
the diagnostics identified in this paper look to the future, pro- 
viding science motivations for next generation hard gamma-ray 
detectors, both of the space-based variety and also ground-based 
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Cerenkov arrays. Since the atmospheric Cerenkov technique is 
intrinsically the more sensitive, due to its much greater collec- 
tion area, the current push to lower their thresholds deep into 
the GLAST band via increased telescope mirror size may in- 
nlucnce the direction gamma-ray burst observational programs 
targeting energies above 1 GeV will take. This advantage is off- 
set by their generally small fields of view (except for water tank 
Cerenkov systems such as MILAGRO), so that large fields of 
view like those afforded by the LAT and the Gamma-Ray Burst 
Monitor (GBM) on GLAST are desirable. 

5. CONCLUSION 

This paper outlines the interesting potential the GLAST mis- 
sion and future hard gamma-ray detectors can offer for probing 
the relativistic nature of outflow in gamma-ray burst prompt 
emission regions. The diagnostics highlighted here use a sin- 
gle, straightforward piece of physics, namely that of electron- 
positron pair creation, and are only dependent on model as- 
sumptions to second order. Spectral attenuation by this pro- 
cess, while not evident in EGRET data for a handful of bright, 
hard bursts, should emerge from the GLAST database, unless 
the bulk Lorentz factors in the burst emission zones are higher 
than presently argued, i.e. greater than around L ~ 10 3 . Then, 
attenuation by interaction of source photons with those of the 
intervening cosmic background field will produce attenuation 
turnovers at above 30-50 GeV, whose energies are independent 
of the GRB flux at a given time. In such cases, bursts can 
be used to explore the evolution of the background field with 
redshift by detectors with sufficient sensitivity in this band. 

For the equally interesting possibility of lower F , pair cre- 
ation turnovers spawned by photons internal to the burst, will 
become prominent in the GLAST energy window. In such 
cases, the representative calculations in Section [31 clearly illus- 
trate how the LAT instrument on GLAST might potentially dis- 
criminate between evolutionary scenarios with constant T , and 
those where F declines with time, should Nature be agreeable 
and offer bright, hard bursts with turnovers below 300 MeV-1 



GeV. A more sensitive, next-generation hard gamma-ray exper- 
iment might additionally be capable of distinguishing between 
adiabatic and radiative evolutionary scenarios: due to their dif- 
ferent r(i) dependence, their turnover energy/flux loci trace 
somewhat different tracks over time. However, this probe will 
be more reliable for flat spectrum bursts and require relatively 
sharp turnovers, or correlated information from the MeV-band 
break evolution. While the GBM on GLAST, in conjunction 
with the LAT, while accumulate a host of data on the MeV- 
band break evolution, next-generation telescopes most likely will 
be needed to differentiate the adiabatic and radiative turnover 
tracks. Should distinction between these two keystone classes of 
blast wave deceleration be possible, the pair creation probes of- 
fer constraints on the explosion energy generated by the central 
engine. Note that while the results presented focus on simple 
broken power-law source spectra, they are readily extended to 
apply to multiple components. They encompass a variety of 
possibilities, and provide a template for burst pair production 
attenuation studies in the soon-to-be realized GLAST era, and 
for future initiatives down the line. 
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APPENDIX 

INTEGRATION OF THE ANGULAR CONTRIBUTIONS TO r 77 (e) FOR STRONG BEAMING 

Here the analytic reduction of the triple integration in Eq. (|15|l is outlined, together with relevant properties of the emergent 
hypergeometric function. It is expedient to scale the integration variables via x u = wx £ and k = x £ p , so that x £ , w and p 
constitute the new set of integration variables. Then, one performs the w integration first, which is easily tractable, and generates 
a simple evaluation in terms of inverse trigonometric functions using identity 2.261 of Gradshteyn & Ryzhik (1980). The x e 
integration is then trivial, absorbing the constraints imposed by the step function <d[Kr)8 m /(2x) ] ■ The result can then be expressed 
(for < f < 1 ) as 

T a « C&W > SM = ^al 2) (jf dy yl+2 ° aTCC0Sy + * 2( ° +2) Si 7 aTCC0Sy } ' (A1) 
with now the one parameter 

* - ~k (A2) 

defining the general character of T a . From this form, the substitution y — sin/? and the invocation of identities 3.621 and 8.335.1 
of Gradshteyn & Ryzhik (1980) lead to the result in Eq. (|16|) when W = 1 . On the other hand, the VP <C 1 limit is readily discerned, 
thereby reproducing Eq. i|18|) ■ These two cases correspond to the results 

G a (1) = A(a) , and ~ * 2 ( 1+ «> , * « 1 . (A3) 
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For '5 > 1 , corresponding to cases where the angular integrations are completely separated form the threshold constraint, Q a (^/) — 
G a {l) — A(a) . For general < ^> < 1 , the integrals for Q a (^f) can be manipulated by the appropriate integration by parts, 
and then the second integral term in Eq. I)A1J| is routinely handled. The first integral is somewhat more complicated, and can be 
developed by using the change of variables y = ^cosx together with the aid of the result 3.671 of Gradshteyn & Ryzhik (1980). 
This renders the first integral expressible in terms of the ordinary hypergeometric function F = 2-F1 , so that 



l+Qt). 



arccos'fr _ gVj - P #F(l/2, a + 3/2; a + 5/2; ^ 2 
a+1 a + 2 + (a + l)(2a + 3) (a + 2) 



< * < 1 



(A4) 



This hypergeometric function reduces to combinations of inverse trigonometric and algebraic functions when a assumes integer 
values, and just algebraic functions when a takes on half integral values. For this range of ^Sf , the hypergeometric function can be 
numerically evaluated using the series expansion 9.100 of Gradshteyn & Ryzhik (1980), yielding 



1 



3 



E 

71 = 



a n z 



2a + 3 
2a + 2n + 3 



n 

fc=i 



2k 



(A5) 



For < z < 1 , this series converges rapidly, faster than a geometric series. In the boundary case of z = 1 , convergence is only 
moderately rapid. The asymptotic behavior of the product in Eq. (|A5(I as n — > 00 can be established by taking its logarithm and 
then using the series identity 44.9.1 of Hansen (1975): 



E 

fe=i 



f- l0 ^( 1 + f) 



-y E x + log e r(x + 1) 



(A6) 



where T(z) is the Gamma function, and 7 E = —ip(l) ~ 0.57721 is the Euler-Mascheroni constant, with ip(z) = d/dz{log e T(z)} . 
By truncating the resultant infinite series and then invoking results 8.365.3 and 8.365.5 of Gradshteyn & Ryzhik (1980), the n — > 00 
limit can eventually be taken, thereby deriving the result 



lim • 

n — >oo 



n / \ 



k=l 



(A7) 



Hence, the series expansion in Eq. (|A5(I converges as n~ 3 / 2 for large n . Such a rate of convergence is not compelling, so that for 
values z ~ 1 it is more expedient to perform the transformation indicated in 9.131.2 of Gradshteyn & Ryzhik (1980), and then use 
the resultant identity 



F 



1 3 5 

2> a + 2'' a + 2' Z 



r(q + 5/2) 
r(a + 2) 



F 



a 



3 I 

2' 2 : 



1 



(2a + 3) Vl^ F 



1, a + 2 



1 



(A8) 



to evaluate the left hand side for domains of (1 — z) <C 1 . The rate of convergence of the two hypergeometric functions on the 
right hand side of Eq. I|A8() is then at least as rapid as a geometric series. This transformation also expeditiously establishes the 
correspondence Q a {^) — > A(a) as $ — > 1~ , using the doubling formula (8.335.1 of Gradshteyn & Ryzhik 1980) for the Gamma 
function. For numerical purposes, it is obviously convenient to use Eq. I|A5(1 for < z < 1/2 cases, and the power series expansion 
about z = 1 of the right hand side of Eq. (|A8|) when 1/2 < z < 1 . 

To conclude, two integrals that are required for the determination of the asymptotic approximation for T(6 mi 77) when r]6 m 3> 1 
are evaluated. Using the integral form for Gai^) in Eq. IjAljl . it can be quickly demonstrated that 



f 1 n ,r\ 



_ 1 - A(a) 
2a 



1 - -4(a) _ 41og„ 2 + l 



4a^ 



8a 



(A9) 



by appropriate integration by parts to generate the derivative dQ a /(M> in each case. 
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